Possibility of Precise Measurement of the Cosmological Power 
Spectrum With a Dedicated 21cm Survey After Reionization 
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Measurements of the 21cm line emission by residual cosmic hydrogen after reionization can be 
used to trace the power spectrum of density perturbations through a significant fraction of the 
observable volume of the Universe. We show that a dedicated 21cm observatory coule probe a 
number of independent modes that is two orders of magnitude larger than currently available, and 
enable a cosmic-variance limited detection of the signature of a neutrino mass ~ 0.05 eV. The 
evolution of the linear growth factor with redshift could also constrain exotic theories of gravity or 
dark energy to an unprecedented precision. 



PACS numbers: 98.80-k, 95.30Dr, 95.55Jz 

Recently, there has been much interest in the feasibil- 
ity of mapping the three-dimensional (3D) distribution of 
cosmic hydrogen through its resonant spin-flip transition 
at a rest-frame wavelength of 21cm[l], Qj- Several experi- 
ments are currently being constructed (such as MWA[33|, 
LOFAR[33,], PAPER 21CMA[3|}and more ambi- 
tious designs are being planned (SKA[3y]) to detect the 
theoretically-predicted emission signal. 

Measurements of the power-spectrum of 21cm bright- 
ness fluctuations could constrain the initial conditions 
from inflation as well as the nature of the dark matter 
and dark energy. The 21cm fluctuations are expected 
to simply trace the primordial power-spectrum of mat- 
ter density perturbations either before the first galaxies 
had formed (at redshifts z > 20) [1, |3| or after reion- 
ization ^ z < Q) - when only dense pockets of self- 
shielded hydrogen (such as damped Lya systems) sur- 
vive 0|. During the epoch of reionization, the fluc- 
tuations are mainly shaped by the topology of ionized 
regions 0, H 0], and thus depend on astrophysical de- 
tails. However, even during this epoch, the line-of-sight 
anisotropy of the 21cm power spectrum due to peculiar 
velocities, can in principle be used to separate the impli- 
cations for fundamental physics from the unknown details 
of the astrophysics 0j[l3|- In what follows, we will focus 
our discussion on the post-reionization epoch [1, [6| which 
offers two advantages. First, it is least contaminated by 
the Galactic synchrotron foreground (whose brightness 
temperature scales with the redshift under consideration 
as (1 -I- z)^'^ [H). Second, because the UV radiation field 
is nearly uniform after reionization, it should not imprint 
any large-scale features on the 21cm power spectrum that 
would mimic cosmological signatures. On large spatial 
scales the 21cm sources are expected to have a linear bias 
analogous to that inferred from galaxy redshift surveys. 
Since a 21cm survey maps the global hydrogen distribu- 
tion without resolving individual galaxies, the 21cm bias 
is expected to be modest compared to surveys that select 
for the brightest galaxies at the same redshifts. 



In general, cosmological surveys are able to measure 
the power-spectrum of primordial density fluctuations, 
P(fc), to a precision that is ultimately limited by cos- 
mic variance, namely the number of independent Fourier 
modes that fit within the survey volume. 21cm observa- 
tions are advantageous relative to existing data sets be- 
cause they access a 3D volume instead of the 2D surface 
probed by the cosmic microwave background (CMB), and 
they extend to a sufficiently high redshift (well beyond 
the horizon of galaxy redshift surveys [ll[) where most 
of the comoving volume of the observable Universe re- 
sides. At these high redshifts, small-scale modes are still 
in the perturbative (linear growth) regime where their 
analysis is straightforward. The expected 21cm power 
extends down to the pressure-dominated (Jeans) scale of 
the cosmic gas which is orders of magnitude smaller than 
the comoving scale at which the CMB anisotropics are 
damped by photon diffusion [3]. 

Altogether, the above factors make 21cm surveys an 
ideal cosmological probe of fundamental physics [1^ . To 
illustrate this point, we show in this Letter that a ded- 
icated 21cm observatory would enable a determination 
of the matter power-spectrum at redshifts z < 6 to an 
unprecedented precision. In our numerical examples, we 
adopt the standard set of cosmological parameters 13 1. 

Number of Modes. The limitation of existing redshift 
surveys of galaxies [3, is apparent in Fig. [1] which 
plots the comoving volume of the Universe out to a red- 
shift z as a function of z. State-of-the-art galaxy redshift 
surveys, such as the spectroscopic sample of luminous red 
galaxy (LRGs) in the Sloan Digital Sky Survey (SDSS) 
[Ul, extend only out to z ~ 0.5 (over ~ 10% of the sky) 
and probe ~ 0.01% of the observable Universe. 

The CMB fluctuations probe a thin shell on the 2D sur- 
face of the sky. The number of modes with a comoving 
wave number k = 27r/A between k and k + dk that fit on 
this 2D surface is, diVcMB — Trkdk [^/(27r)^] , where A = 
D'^dVl, D is the comoving distance to the last scattering 
surface at z ~ 10^ and dfi is the solid angle of the sur- 
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FIG. 1: The fraction of the total comoving volume of the 
observable Universe that is available up to a redshift z. 



vey field. Redshift surveys of galaxies and 21cm surveys 
probe a 3D comoving volume V and potentially access a 
larger number of modes, dN-^^^ — 2'Kk^dk [V/(27r)'^] . 

Figure [5] compares A^cmb = {k / 10) dNcMB / dk with 
-^3D = {k / 10) dN3Y)/dk, for future 21cm surveys after 
reionization [5|. The CMB data set is assumed to cover 
a fraction /sky = 0.65 of the sky (excluding the re- 
gion around Milky Way galaxy). For comparison, we 
also show the corresponding number of modes within the 
same k interval in the spectroscopic LRG sample of SDSS 
ll|, which covers ^ 3700 square degrees out to z ~ 0.5 
or a volume of V = 0.72/i^'^Gpc'^ (where h « 0.7 is the 
Hubble constant in units of 100km s~^ Mpc~^). 

21cm observatories that are currently under construc- 
tion (such as MWA) will survey only a few percent of the 
sky and process only ^ 15% of the available frequency 
range (band-pass). In this Letter, we consider future 
21cm surveys that would potentially cover /sky — 0.65 
with a processed frequency-bandwidth spanning a red- 
shift range of a factor of 3 in (1 -I- z) centered on z = 1.5, 
z = 3.5 and z — 6.5 With an array design based on 
MWA in which the effective area of each tile of 16 dipole 
antennae equals its physical area, the value of /sky = 0.65 
corresponds to 33 correlated primary beams or fields. 

Results. The fractional uncertainties in the 21cm 
power-spectrum P21 for a cosmic- variance limited survey, 
{AP21/P21) = 1/VN, are presented in the inset of Fig. [3] 
(straight lines). Also shown in both the inset and main 
panel are the noise curves for observations using a design 
based on the so-called MWA5000 experiment, which is 
assumed to have 10 times the collecting area of MWA 
(as described in Refs. 0,lii). MWA5000 would be cos- 
mic variance limited in an integration time of ~ 10'^ hours 
at wave numbers near k ~ O.lMpc"^. Since the dipoles 
of each antenna tile look at ~ tt steradians, simultane- 
ous processing of multiple primary beams would allow a 
survey of area /sky = 0.65 at 10^ hours of integration per 
pointing within a few years. In computing the thermal 
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FIG. 2: The number of modes A'^ within a wave number bin of 
width Afc = fc/10 centered on fc, that are available in different 
cosmological surveys. The thick dashed grey line corresponds 
to the spectroscopic LRG sample of SDSS '11'], while the thick 
solid line (marked HRG) corresponds to a future spectroscopic 
survey at 2.5 < z < 3.5 covering 1000 square degrees with a 
a co-moving galaxy density equal to the LRG sample. The 
thick dark line corresponds to a CMB data set with /sky = 
0.65. The thin lines show the number of modes accessible 
in a 21cm survey (including the limit on large scale modes 
due to foreground removal Q) covering /sky = 0.65 within a 
redshift range spanning a factor of 3 in (1 -)- z), and centered 
on z — 1.5,3.5 and 6.5. For z < 1.5, we have truncated the 
curves at k = 0.2 Mpc~^ to illustrate the smaller range of k 
accessible within the linear regime at lower redshifts [2a |. 



noise of the observatory, we have adopted a model for the 
bias factor &21 of the 21cm sources (63]^ = P21/P) from 
Ref. The noise curves include a limit on large scale 
modes due to foreground removal. Current estimates for 
MWA show that foreground removal should be effective 
for modes over a frequency range < 6MIIz which is ^ 
of the total 24MHz processed bandwidth However, 
improvements on this range would provide access to a 
larger N as well as to lower-A; modes. In Fig. [3] and sub- 
sequently, we assume a scenario in which the foreground 
can be removed on scales of up to of ^^'^ total processed 
bandwidth, namely x (V3-l/vl)[1400MHz/(l-hz)] = 
30[(H-z)/4.5]-iMHz. 

For comparison, we also show the noise curves for 
the SDSS-LKG survey (thick grey line), including the 
effects of Poisson shot-noisc 38] , (APgai/Pgai) = [1 -I- 
(&^P(fc)rigai)^^]/\/]V, for a galaxy number density of 
rigai = 46748/(0. 72/i~3)Gpc"^ and a bias factor of 6gai = 
2 11|. We find that the potential 21cm constraints on 
the matter power-spectrum are 1-2 orders of magnitude 
better than a low-redshift galaxy survey like SDSS-hKG. 

As an example for the potential use of a 21cm sur- 
vey, we show in Fig. [3] the expected relative changes 
in the amplitude of the power-spectrum owing to the 
presence of a massive neutrino [l8|. At wave numbers 
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FIG. 3: The fractional change in the amphtude of the power- 
spectrum owing to the presence of a massive neutrino (hor- 
izontal grey lines, asymptoting towards a constant at high 
k values). The case shown, = 0.004, corresponds to 
— 0.05eV. For comparison, the limits imposed by cos- 
mic variance on measurements of the power-spectrum from 
the SDSS-LUG and a future 1000 square degree galaxy sur- 
vey at2.5<2:<3.5 are marked by the thick dashed and solid 
grey lines respectively. The [/-shaped error curves correspond 
to an all-sky 21cm survey [/sky ~ 0.65 over a redshift range 
spanning a factor of 3 in (1 -|- z)] with MWA5000 and a 10^ 
hour integration per field (line styles for z — 1.5, 3.5, 6.5 as in 
Figure [2}. The noise is plotted in bins of size Ak/k = 0.1. 
The inset shows these results on logarithmic axes that span 
a larger dynamic range of achievable precision. The straight 
thin lines in the inset show the cosmic-variance uncertainty 
in the power-spectrum measurement owing only to the num- 
ber of available modes. For z < 1.5, we have truncated the 
curves at = 0.2 Mpc~^ to illustrate the smaller range of k 
accessible within the linear regime at lower redshifts [2a |. 



much larger than the neutrino free-streaming wave num- 
ber {kfs = 0.lilm{0)hy/J^ Mpc~-M, the suppression of the 
power-spectrum is given by 2^, 21 1, 



PjkJ.) 
P(fc,/. = 0) 



1.9 X 1055(0)17™(0) 



A. 



(1) 

where g{z) « + 0.003(r2A/f7™)'*/^] is the growth 

function of the gravitational potential for matter and vac- 
uum density parameters of U,m{z) = [1 + f^A/i^m(0)(l + 
and f^A = (1 - ^m) 0; fu = = 
0.087Vi,(my/l eV) is the present-day mass fraction of the 
matter density carried by Ni, neutrino species of particle 
mass rriu. In this Letter we conservatively assume a non- 
degenerate hierarchy of neutrino masses with Ni, = 1 and 
rriu denoting the largest mass eigenstate. Figure [3] shows 
the case corresponding to fi, — 0.004 {niy = 0.05eV). 

While Fig. [3] indicates that the cosmic variance in a 
galaxy redshift survey (such as the SDSS-IMG survey) 
is sufficiently small to detect the suppression of power 
at fc ^ kfs due to a neutrino mass of mv = 0.05eV, un- 



certainties in other cosmological parameters and param- 
eter degeneracies reduce this sensitivity by an order of 
magnitude [isl. fioj . Thus, in order to avoid possible sys- 
tematic offsets between the power-spectrum amplitude 
observed by different techniques (such as galaxy surveys, 
CMB maps, and Lya forest data) at different k values, 
it is desirable for the 21cm survey to be self-contained 
and cover a sufficiently broad range of wave numbers 
that probe the curvature of the neutrino effect in Fig. 
[3l For this to be achieved, the removal of the Galactic 
synchrotron foreground would need to be effective over 
large frequency intervals of up to ~ 30[(l-|-2;)/4.5]~^MHz. 
Foreground removal across such intervals would provide 
access to the required range of scales over which a con- 
stant offset (in the form of a linear bias factor) would 
not affect the rrii, measurement. The desired wave num- 
ber to be reached by foreground removal corresponds to 
the scale where cosmic variance is larger than the change 
induced by a massive neutrino. Currently, the detailed 
properties of the foreground are not well measured. The 
feasibility of foreground removal over a broad frequency 
interval will remain uncertain until data from the first 
generation of 21cm observatories is analysed. 

Figure [3] illustrates that foreground removal to ~ 
30[(1 -f z)/4.5]"iMIIz would be sufficient to detect the 
modification of the power spectrum due to the minimum 
neutrino mass of 0.047 ± O.OleV indicated by the latest 
atmospheric neutrino data [1^, [l^] at all spatial scales 
where the effect is larger than cosmic variance. The ad- 
vantage of the large survey volume is evident since it 
allows the modification of shape to be measured in addi- 
tion to the suppression detected by the SDSS-IMG sur- 
vey (the latter being cosmic variance limited on scales 
where the shape is measured). In a follow-up paper 
[2^ . we will address the level of degeneracy with other 
cosmological parameters. Already, Fisher-matrix stud- 
ies 0, [l^l have demonstrated the improved capabilies 
of 21cm observations during the epoch of reionization, 
where contamination from astrophysical sources needs to 
be removed through the angular dependence of the 21cm 
power-spectrum 

A 21cm survey measures the modulation in the cumu- 
lative 21cm emission from a large number of galaxies, 
as its coarse angular and redshift resolution is not ca- 
pable of resolving the 21cm sources individually (Bl. [l^. 
The damped Lya systems which contain most of the hy- 
drogen mass in the Universe at 2; < 6, are expected to 
be hosted by abundant low mass galaxies [2^ and thus 
have a weak bias relative to the underlying matter dis- 
tribution on large spatial scales. This weak bias is not 
expected to introduce a feature to the power-spectrum 
that is degenerate with the neutrino signature (as would 
be the case prior to reionization). Although the comov- 
ing wave number at which non-linear evolution becomes 
important increases from k ~ 0.1/i~^Mpc~^ at z = 0.3 
to ^ Q.bh~^Mjpc~^ at z = 3 [2^, the constraints on rrii, 
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can be potentially improved by accounting for the related 
non-linear effects [28{ . 

By measuring the evolution of the growth factor with 
redshift to the exquisite precision implied by Fig. [31 it 
would also be possible to constrain alternative theories 
of gravity or dark e nerg y well beyond the capabilities 
of existing data sets [29|. The evolution of the growth 
factor would be limited by the knowledge of the bias 
of the 21cm sources, 621. The limit on the uncertainty 
in growth factor would satisfy [d\ng/d\n(l + z)] « 
[dln62i/dln(l + z)]. 

Finally, we note that a precise P{k) measurement 
at multiple redshifts would also allow to determine the 
redshift evolution of the baryonic acoustic oscillations 
(BAO) in the 21cm power spectrum I s l . 1301 . The BAO 
scale constitutes a standard ruler [26|, |27| that can be 
used to measure the equation of state of the dark energy 
S [3 ) constrain 1 — i7tot , and further remove degenera- 
cies between mi, and other cosmological parameters Isij . 

Hardware. The required 21cm observatory could be 
similar in antenna design to the planned MWA but would 
require expanded collecting area and computer resources 
to account for the increased cross-correlation require- 
ments and the analysis of multiple beams. A suitable 
observatory would contain ten times the number of tiles 
in MWA. The computational load increases as the square 
of the number of tiles in a telescope, and linearly with 
both the amount of processed band-pass and the num- 
ber of cross-correlated primary beams. Thus, an all-sky 
survey over a frequency range covering a factor of 3 in 
requires an overall improvement by ^ 10^ in com- 
puter power relative to MWAjsoj. 
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The factor of 3 in (1 -I- z) corresponds to the largest fre- 
quency bandwidth over which a low-frequency dipole an- 
tenna has suitable sensitivity. 

We ignore Poisson fluctuations in the 21cm power spec- 
tra since the contributing galaxies are expected to be of 
much lower mass than LRGs, and there should be a large 
number of them per resolution element of the survey. 
The cost of the antenna hardware is expected to be in the 
range of hundreds of million of dollars, an order of mag- 
nitude higher than for MWA. Improved telescope designs 
[2J and advances in computer technology could help to 
realize our projected 21cm observatory. 



